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• Abstract 

Based on extraction of resonances from quantal time delay, a theo- 
^\ ' rem relating quantal time-delay and the number of resonances below a 

04 ■ certain energy is proved here. Several illustrations from quantum me- 

chanics, neutron reflectometry and hadron resonances are presented. 
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CN \ Physics of weakly-bound systems and resonances has been of great inter- 

est for the important role it plays in nuclear and particle physics. In particu- 
(T) . lar, deriving reliable information about unstable and short-lived states leads 

to a deeper understanding in the theory of elementary particles in case of 
T ^ ' hadron resonances and hypernuclei [2], compound nuclear resonances jH] 

in (n, 7) and (p, 7) reactions and in development of new models and theories 
in the active field of radioactive nuclei [I] . The research on halo nuclei is also 
intimately connected with the search for superheavy nuclei |S]. The struc- 
^ | ture of these elusive states is usually in the form of two-body, three-body or 

many-body resonances It would be interesting if we could infer from the 
. I phase-shift data for ir-ir, or, r]-d, or, a-a scattering whether respectively a 

^ ■ p-meson around E = 770MeV is formed, or, a hyper-nucleus, or, a resonant 

a-cluster (7j is formed. Clearly, we would like to be sure of our analysis by 
studying simple, non-trivial illustrations from quantum mechanics. In this 
paper, we present various illustrations to elucidate how resonant structures 
can be reliably extracted and understood. 

The concept of central importance for our considerations is that of time- 
delay. Time-delay, T is defined as the difference between the density of 
states with and without interaction. This reduces to the following relation 
[5] between T and the S'-matrix [S], 

T = -^trSt-^| (1) 
dE 
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for a general form of interaction. For central potentials, for the l-th partial 
wave, Ti reduces to fidb\jdE where <5/ is the difference between phase shift 
5i(E) and hard-sphere scattering phase-shift 5j^{a,E), and where a is the 
range of interaction. For an s-wave scattering, 5l I {a : E) = —ka where k is 
the wavenumber. 

The complex-energy poles of time-delay on un-physical sheet like those 
of S'-matrix represent resonances. Consequently, the peaks in T{E) are 
signatures of resonances jH]. 

In a system where number of single-particle states modify in the presence 
of an impurity, the change in the states is related to the energy-derivative 
of phase-shift. The total change in the density of states leads to the Friedel 
sum rule, which is basically a statement of charge-neutrality On the 

other hand, an energy-integral of the phase-shift is shown to be connected to 
the energy of the impurity - Fumi's theorem. All the filled states contribute 
to the energy of the impurity via electron- impurity interactions. 

In the context of nuclear physics, some recent works treat the effect 
of continuum by considering the contribution of naroow resonances. These 
narrow resonances are assumed to be like broadened bound states. This is, 
of course, not a valid assumption for broad resonances occurring in parti- 
cle physics. When this assumption holds, the change in the single-particle 
continuum level density, 5p, in the presence of a potential V is found by 
subtracting the free particle level density po from the total level density, p : 



,„ r| _„„ ) = i; w±i2i!ja, 



the last being due to Friedel ^2j. The right hand side of (2) is connected with 
time-delay. We will consider an energy-integral of time-delay in this paper 
and illustrate it by a variety of physical situations of simple and advanced 
nature. 

For a narrow resonance at energy Eq of width T, time-delay has a 
Lorentzian or a Breit-Wigner (B-W) form : 

T( E) = hT/2 „. (3) 

Then, it follows that the number of resonances below E* for a fixed partial 
wave is given by ^3] 

- - is: m - 

= N + A, < A < 1. (4) 

This is related to the approaches based on level-density E]- Since 
most of the methods fail to determine the absolute values of phase-shifts, 
therefore, R.H.S. of (4) would be different from the expression [5i(E*) — 
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6i(Q)]/ir. However, following somewhat more general considerations along 
the lines of [^j , the spectral property of time-delay implies that 



This is valid also for broad resonances, as shown in the Illustration IV be- 
low. To appreciate (5), following remarks are in order : (i) (5) is not the 
Levinson's theorem; (ii) if the energy-range includes bound states also, sin- 
gularities will appear for bound states and non-zero widths appear for res- 
onances, shapes being Lorentzians only for narrow ones (see also comments 
at the end of Illustration II below). 

Construction of time delay above using the complex energy-eigenvalues 
may remind one of the potential scattering formalism due to Kapur and 
Peierls |16| . To illustrate, with the s-wavefunction in the range of potential, 
the complex roots, £ of hdu/dr(£,r = a) = W2m£u(£ , r) are taken as 
resonances here. However, in |16j . y/2m£ is replaced by y/2mE, where E is 
fixed. Moreover, Kapur- Peierls' formalism is designed to yield a non-unitary 
S'-matrix. 

At resonance, the phase-shift 5(E) = arctan ( ) , due to the pres- 
ence of an inverse trigonometric function, has an essential ambiguity of nir 
where n remains uncertain. On the other hand, the ambiguity disappears 
in T(E) as one gets a purely algebraic expression, making time-delay un- 
ambiguous and hence more reliable than the phase-shifts themselves. The 
second derivative of 5(E) reveals that on the either side of E = Eq the cur- 
vature changes sign and hence a point of inflexion in the variation of 5(E) 
is a necessary signature for the existence of a resonance. 

Resonances which are complex energy Gamow-Siegert states of a poten- 
tial generally arise when the potential is localized (e.g., short-ranged), or, 
when there is a barrier of a finite width attached to the well. The square 
well represents the former when I = and the latter when I > 0. However, 
for simplicity we present s-wave square well potential [T7] and s-wave Dirac 
barrier potential Vq5(t — a) ^Tf\ which represents the second situation. It is 
worth noting that the Illustrations III and IV discussed below are novel as 
they do not fall in the above categorization. Believing that a fundamental 
result is understood better in terms of diverse examples, we present them 
below. 

Illustration I : Square- well potential (for barrier, see |18j ) 

To begin with, we consider the most familiar potential V(r < a) = 
—Vq, V(r > a) = 0, as an example to illustrate our results on resonances. It 
may be noted that a reliable description of narrow resonances is important 




(5) 
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for separating the background. The s-wave S— matrix can be written as 



So(E) = 



ik tan pa + p 



(6) 



ik tan pa — p 



The complex poles (resonances) of So(E), i.e., the roots of ik tan pa = p, 
are equivalent to demanding outgoing wave, uq{t > a) = e tkr , instead of 
uq(t > a) as given above. The s-wave time delay for square- well potential 
can be obtained as 



Relevant complex poles of S-matrix, (Ej — iTj/2), such that Ej,Tj > 
0, are known to be the resonant energies. With the help of FindRoot by 
Mathematica, by taking 2m = 1 = h, and fixing the values of Vo and a, we 
find the resonant energies (e.g., poles of (6), for / = 0). We then construct 
time-delay as a sum of Lorentzians as these resonances are narrow. The 
solid lines in Fig. 1 display the exact values of %{E) for I = 0,1,9,10 
for square-well potential (Vq = 5, a = 10) up to an energy, E = 10. The 
dashed lines show the time-delay as calculated by using first fifteen complex 
energies Ej — iTj/2. Excellent reproduction of the time-delay can be seen 
to testify the representation of the time-delay by a sum over Lorentzians 
for a localized potential possessing resonances. We check that Ej are the 
energies where the peaks in time-delay occur. Further, we check that Tj is 
nothing but two times the inverse of the height of the j th -pe&k. The values 
of the integral of time-delay up to E = 10, ur (4), are mentioned in each 
part of Fig. 1 which can be verified to be giving the number of resonances 
as depicted by the number of peaks in Ti(E). 

Interestingly, in the above calculations when I = 9, 10, we find that the 
first complex roots E = 0.38499 - 0.479894i = E A and E = 0.541725 - 
0.574161? = Eb respectively do not show up as peaks in T\(E). Moreover, 
if they are included in the summation (4), they spoil the reproduction of 
the time-delay at other energies. Similarly, when Vo = —5 and a = 2, we 
find that E x = 0.023387 - 0.542466i = E c and E 2 = 9.38365 - 4.43007i 
are the poles of Sq(E) (6). We find that E\ does not lead to a peak in 
Tq(E), whereas E 2 shows up as a broad peak in Tq(E). Notice that in these 
examples Tj > Ej and the corresponding B-W profile appears only partly 
for energies, E > 0, most of it creeps in the negative energy regime. We rule 
out Ea, Eb, Ec as resonances by studying the corresponding wavefunctions. 
The wavefunctions at energies equal to the real part of these complex en- 
ergies show scattering-state-like behaviour wherein the localization within 
r < a is inhibited. 

Illustration II : Repulsive Dirac delta s-wave barrier 



%{E) = h 



Vq tan pa + apk 2 sec 2 pa 



(7) 



2pk(p 2 + k 2 tan 2 pa) 
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Dirac delta barrier, Vq5{t — a), Vq > in three dimensional Schrodinger 
equation ^7j, even in s-wave, supports the Gamow-Siegert states. The poles 
of S-matrix can be worked out to be the complex roots of k cot ka — aVo = 
ik, a condition which one also gets by using the outgoing wave boundary 
condition at r = a. The s-wave quantal time delay can be found as 

aVo tan 2 ka + ak 2 sec 2 ka 
o{ - ' = 2[k 2 tan 2 ka + (k + aV tan ka) 2 } ' ^ 

For Vq = 10 and a = 1, we find first four poles of S-matrix and construct 
time-delay as a sum of Lorentzians (shown as dotted lines in Fig. 2). Even 
with only four poles, the Tq(E) (8) (solid line) is excellently reproduced by 
the dotted line in Fig. 2. By an integration of time-delay up to energy 
E = 170, we get the value of the integral ur = 4.0114, correctly indicat- 
ing four resonances or metastable states. In the limit of Vq — > oo, we get 

'2 2 ^2 

Ej — ► J 2 ma i anc ^ Tj — ► these states will be the bound states of a particle 
contained between two rigid walls of width a. Thus, for large values of Vo 
one gets Dirac delta like spikes in Tq(E) at energies E = Ej. 

Illustration III : 

While studying the reflectometry of the polarized neutrons from magne- 
tized superconductors, Zhang and Lynn |19j arrived at a novel result wherein 
a single, pronounced, parameter-dependent dip occurs in the reflectivity. In 
terms of simple quantum mechanics this means that the reflectivity, R(E), 
of the potential-step barrier, V(x) = V\ + V^l — e~ x / a ), as a function en- 
ergy shows a single, pronounced, spike- like minimum at an energy slightly 
above the step, for certain sets of parameters, (Vi,V2,a). Usually when a 
semi-infinite potential is smooth, R{E) is a smoothly converging function of 
energy. A discontinuity of derivative at a point in a semi-infinite potential 
turns out to be the essence [201 121) of this novel reflectivity. An interesting 
claim (221 that the energy, (E = E^), at which this reflectivity minimum 
occurs corresponds to a 'half-bound state' has been refuted by the help of 
several model potentials |2j|. Let us take the case of an exponential step 
where r(E) = ' R(E)e if) ( E ' is given by 

r(E) _ ikJ-2ipa(2ga) + qJ'_ 2ipa {2qa) 

n } ~ ikJ_ 2ipa (2qa) - qJ'_ 2ipa (2qay { } 

where p = V 2m ( E ~ Vl V2 ) anc [ ^ _ ^J 2 ™Y 2 . Whenever 2qa is close ^5] to the 
zeros of Jo(z), a single pronounced minimum exists in R{E). More precisely, 
when V\ = V% = 1 and a = 1.31 (2m = 1 = K) the dip in reflectivity occurs 
at E = 2.0445 = E d j5U]. We find that reflection time-delay, possesses 
a single dip at, E = Ed, confirming the existence of a resonance there (see 
Fig. 3). Since the contention is regarding the dip in reflectivity, we have 
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calculated the reflection-phase-shift (6) and its derivative. Eventually, the 
integration (4) of time-delay over a long range of energy (from E = 2 to 
E=4,6,10) divided by tt has been checked to yield a number very close to 1. 
Thus, the composite semi-infinite potential steps |201 121j are a new type of 
model potentials which are neither localized nor there is a well attached to 
a barrier. They, however, possess a single resonance. We also find that the 
low energy, single minimum in reflectivity of the composite potential wells 
[2*0] such as V[x) = — Vq exp(— |rc|/a) is nothing but a resonance. 

Illustration IV : 

Usually, narrow resonances are treated as elementary particles [23j . How- 
ever, even relatively broader resonance like the excited state of a nucleon 
called A ++ (1232) are treated as particles. Here we would like to demon- 
strate that time delay calculated from the available J21j elastic scattering 
phase shifts of tt + and p (P33) yields a peak at E = 1218 MeV of width 
r = 129 MeV (see Fig. 4) corresponding to the formation of the well- 
known baryon resonance called the A-resonance which has mass as 1232 
MeV and width as 120 MeV. This result can also be found in "2H] where 
there is a more detailed study on unflavoured baryon resonances. The inte- 
gration (4) of time-delay in the energy range given in Fig. 4 divided by tt 
yields 0.87 instead of unity, simply because the resonance data is truncated 
at the lower energies. This illustration, in addition to "2""] shows that (4) 
and (5) are also valid for non-central interactions and for broad resonances 
with lifetimes ~ 10~ 23 s. It should also be noted that one is in a relativistic 
regime with energies of the order of GeV. The concept of time-delay holds 
in relativistic regime because of its connection with the 5-matrix (1) [26j . 

Hitherto, the energy-integral of the first derivative of phase shift over a 
long range of energy is supposed to either vanish [37] or yield the number 
of bound states [2B] via the well-known result in (10) below. In such a 
confusing scenario, we have brought out the correct meaning of the energy- 
integral mentioned above. We have revealed that the integral (4) does not 
vanish and it yields the number of resonances possessed by a potential. 
Since time and energy are conjugate variables, (5) can be interpreted as 
an analogue of the Bohr-Sommerfeld quantization. This is also a quantum 
analogue of some results in classical scattering theory obtained twenty years 
ago by Narnhofer and Thirring 28 . 

Finally, we would like to make some observations on the Levinson's the- 
orem in its original form relates the phase-shift at zero energy to the number 
of bound states, ne, possessed by an attractive potential : ns = 5i(0)/tt 
[HI EH]- However, it is very often written as [SI "231 EDJ ED] 

n B = MO) - <5/(oo)]/7r. (10) 

It is important to note that (10) seems to relate the negative-energy bound 
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states to the phase-shift at an infinitely positive energy. This is basically 
done to provide a reference at infinite energy where phase shift is assumed 
to be zero. Although for short-range potentials, <5/(oo) is zero, this is not 
generally true [3U|. In a specific calculation, one may find <5;(0) and <5;(oo) 
with ambiguous factors of, say, rayix and rri2'K with mi,m.2 as two arbitrary, 
uncontrollable integers, having nothing to do with the potential. To obtain 
phase shifts free of such ambiguities (modulo 7r), one needs to employ special 
methods like the variable phase approach [221 an d an integral representation 
of phase shift found recently [33] • 
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Figure Captions : 

Fig. 1(a): Time-delay, T(E), for a square well-potential {Vq = —5, a = 10) 
for four values of angular momentum / = 0,1,9,10. The solid lines are 
the exact calculation and the dashed lines represent a sum of Lorentzians 
using first fifteen resonances (poles of S'-matrix on un-physical sheet). In- 
terestingly, in the case of I = 9, 10 the very first complex energies Ea = 
0.38499 - 0.47984z and E B = 0.541725 - 0.57416H, respectively are not in- 
cluded (see text for details). Notice that the values of tir agree with the 
number of peaks in all the plots of time-delay, illustrating (4) . 

Fig. 1(b): The first peaks corresponding to Fig. 1(a) are separated out 
and shown for clarity. 

Fig. 2: Excellent reproduction of the time-delay (8) by a sum of Lorentzians 
for the Dirac delta potential by using just first four complex energies (res- 
onances). A value of tir close to 4, illustrates existence of four resonances 
up toE = 170. 

Fig. 3: An instance of reflectivity, R(E), dip in the exponential semi- 
infinite potential step when V± = V2 = 1 and a = 1.31. At E = 2.0445 
there occurs a dip in R(E) and T(E) and the (reflection) phase-shift, 9(E) 
displays point of inflexion. The value of tir from £ , = 2upto.E' = 4,6,10 
has been found to be very close to 1. 

Fig. 4: The well known baryon resonance, A ++ is demonstrated as a 
peak in time-delay in the phase-shift data |24| of the scattering : ir + + p — > 
A++. We get M A = 1218 MeV and T A = 129 MeV, whereas the standard 
values are M A = 1232 MeV and T A = 120 MeV. The value of n R we get is 
0.87. 
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